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Abstract 

We have performed detailed multifractal analysis on the minutely volatility of two indexes 
and 1139 stocks in the Chinese stock markets based on the partition function approach. 
The partition function Xqi s ) scales as a power law with respect to box size s. The scaling 
exponents r(q) form a nonlinear function of q. Statistical tests based on bootstrapping 
show that the extracted multifractal nature is significant at the 1% significance level. The 
individual securities can be well modeled by the p-model in turbulence with p = 0.40 ± 
0.02. Based on the idea of ensemble averaging (including quenched and annealed average), 
we treat each stock exchange as a whole and confirm the existence of multifractal nature in 
the Chinese stock markets. 
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1 Introduction 



Since the pioneering work in the 1990's [Q], Q, yfl, there has been a vigorous con- 
tinuing investigation aimed at discovering remarkable similarities between finan- 
cial markets and turbulent flows. Multifractal analysis, which was initially intro- 
duced to investigate the intermittent nature of turbulence has also been 
extensively applied to various financial time series [7]. Multifractality has been 
regarded as one of the most important stylized facts in equity returns. Many dif- 
ferent methods have been applied to characterize the hidden multifractal behav- 
ior in finance, for instance, the fluctuation scaling analysis fiji, 101, the struc- 



IJ, lis O, LL8L [JJ, |20j, |21|], multiplier methoaj 22JL multifractal detrended flue 



tuation ana 



041l,|42j,|43j 



ysis (MF-DFA) |23, 24, 



method J3lL |32j, |3j, |34|, |35, 



36, 



37 



25 
38 



26, 
39, 



ture function (or height-height correlation function) method all! 



27, |28l |29|, |30J], partition function 
4011 . wavelet transform approaches 



44J, to list a few. 



To the best of our knowledge, the first application of partition function approach 
was carried out by Sun et al. to analyze the multifractal nature of minutely Hang 
Seng Index (HSI) data for individual trading days from January 3, 1994 to May 28, 
1997 113 lh . Ho et al. also reported that there exists multifractal feature in the daily 
Taiwan Stock Price Index (TSPI) data from 1987 to 2002 H- Wei and Huang an- 
alyzed the 5-min intraday data of the Shanghai Stock Exchange Composite (SSEC) 
index for individual trading days from January 1999 to July 2001 [34]. They found 
a different empirical result from the HSI and constructed a new measure of mar- 
ket risk to predict index fluctuations. A similar multifractal feature was reported in 
the high-frequency data of the SSEC index at different timescales and time peri- 
ods 11361 13711 . What is more interesting is the claim that the observed multifractal 
singularity spectrum has predictive power for price fluctuations and can serve as a 
measure of market risk I132l l34l. 13811. 



The partition function approach was introduced to characterize singular measures 
[|45|], where a spectrum of local singularities exist. When the measure is homoge- 
nous in singularity (thus not a multifractal), the scaling exponent r(q) of the g-order 
partition function is linear with respect to q. In all the aforementioned studies, the 
index per se is treated as being proportional to a singular measure. Intuitively, there 
is no singularity in the index. This argument is indeed confirmed by the very nar- 
row width of the extracted singularity spectrum f(a). More precisely, the minimum 
singularity strength a min and the maximum singularity strength a max are both close 
to 1. Alternatively, the scaling exponent function r(q) is linear against q, which dis- 
approves the presence of multifractality. Zhou discussed these issues using 5-min 
SSEC data and performed extensive statistical tests [391]. Jiang and Zhou further 
clarified the situation using intraday high-frequency data for HSI, SZSC (Shenzhen 
Stock Exchange Composite), S&P500, and NASDAQ H. 
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Having said this, we nevertheless figure that the idea is valuable to apply parti- 
tion function approach to the multifractal analysis of stocks and indexes and to the 
possible application of multifractal properties in market prediction and risk man- 
agement. This, of course, should be done using returns rather than stock prices 
or indexes. In this work, we shall follow this line to perform detailed multifractal 
analysis on Chinese stocks and indexes, which form a huge database. In almost all 
previous studies concerning the multifractal nature of financial markets, the investi- 
gation was conducted based on a single financial time series for individual stocks or 
indexes. In this paper, we propose to study many stocks as an ensemble. The idea is 
that different stocks in a same market share many common underlying mechanics. 
This consideration leads to the assumption that different stocks are realizations of a 
same stochastic process and we can perform ensemble av eraging in the multifractal 



analysis as an analogue of diffusion-limited aggregation 11461 1471 l48l. 14911 



This paper is organized as follows. In Section [2} we describe the data sets and 
perform preprocesses. The methods adopted are explained in detail in Section [3] 
Section HJpresents the results and discussions. Finally, Section [5] concludes. 



2 Preliminary information 



2. 1 The date sets 



We use a nice high-frequency database in the Chinese A-share stock markets from 
January 2004 to June 2006. The trading rules were not changed during this time 
period. After eliminating those stocks that have recording errors or less than 0.1 
million data points, we are left with 715 stocks listed on Shanghai Stock Exchange 
(SHSE) and 424 stocks traded on the Shenzhen Stock Exchange (SZSE) (1 139 in 
total). We also include two indexes, the Shanghai Stock Exchange Composite Index 
(SSEC) and the Shenzhen Stock Exchange Composite Index (SZSC). The average 
size for each stock is about 1.32 million and there are totally more than 1.5 billion 
data points for all the stocks. 

All of the these indexes and stock prices are tick-by-tick data, which were recorded 
based on the market quotes disposed to all traders in every six to eight seconds. 
Each datum is time stamped to the nearest second at which one transaction occurs. 
Hence, the recording time interval of price series is uneven for each stock. For 
convenience, denote the time sequence for each data set as i , and the corresponding 
price sequence as p(ti), where i = 1, 2, • ■ ■ . In literature, i is also called the event 
time. 
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2.2 Definition of volatility 



For each time series, we define the event-time return r(tj) over one event step as 
follows 

r(t i ) = \np(t i )-\np(t i _ 1 ) . (1) 
We then calculate the minutely volatilities as follows, 

V (t)= £ [r(r)l , (2) 

t-At<T^t 

where At = 1 min. On each trading day, the Chinese stock markets contain opening 
call auction (9:15 a.m. to 9:25 a.m.), cooling period (9:25 a.m. to 9:30 a.m.), and 
continuous double auction (9:30 a.m. to 11:30 a.m. and 13:00 p.m. to 15:00 p.m.). 
Since July 1, 2006, the Shenzhen Stock Exchange introduced closing call auction 
(14:57 to 15:000. Our analysis focuses on the data recorded during the continuous 
double auction. In this way, there are 240 data points of the volatility for each time 
series on a single trading day. 



3 Methodologies 



3.1 Partition function approach 



Denote the minutely volatility series as {v(t) : t = 1, 2, • • •, T}. Then the series is 
covered by N boxes with equal size s, where N = [T/s]. On each box, we define 
a quantity u as follows, 



u\n 



;s) = u([(n-l)s + l,ns}) = - l)s + t] , (3) 



t=i 



where \{n — l)s + l,ns] is the n-th box. The box sizes s are chosen such that 
N = [T/s] = T/s. The measure fj, on each box is constructed as follows, 



u(n; s) 
£^ =1 w(m;s) 



/x (i n, o> 

M^) = , , (4) 



We then calculate the partition function x q 1)4511 



x q (s) = EK«)]' > ( 5 ) 



n=l 



1 During the time period of our data sets, both exchanges did not have closing call auction. 
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and expect it to scale as 

Xq (s) ~ s t(<?) , (6) 

where the exponent r(q) is a scaling exponent function. The local singularity ex- 
ponent a of the measure y, and its spectrum f(a) are related to r(q) through a 
Legendre transformation Il45ll 



a = dr(q)/dq 
/(«) = 9« - r (<?) 



(7) 



When w(n; s)/Yl u(m; s) C 1 and q 3> 1, the estimation of the partition function 
X will be very difficult since the value is so small that the computer is "out of the 
memory." To overcome this problem, we can calculate the logarithm of the partition 
function, ln% g (s), rather than the partition function itself. A simple manipulation 
of Eqs. © and © results in the following formula, 



In Xg(s) 



N 

71=1 



u(n; s) 



max{w(m; s)} 



+ gin 



niax{«(m; s)} 



(8) 



where max{w(m; s)} is the maximum of u(m; s) for m — 1, 2, 



N. 



3.2 Bootstrapping for statistical test 



To test for the possibility that the empirical multifractality could be artifactual, we 



adopt the following bootstrapping approach 11391 14011 . For a given volatility time 
series, we reshuffle the series to remove any potential temporal correlation and 
carry out the same multifractal analysis as for the original data. We impose a very 
strict null hypothesis to investigate whether the singularity spectrum f(a) is wider 
than those produced by chance. The null hypothesis is the following: 

H l : Aa ^ Aa rnd . (9) 

The associated probability of false alarm for multifractality (so-called "false posi- 
tive" or error of type II) is defined by 

= #[A a <A a „ J ] i 
n 

where n is the number of shuffling and #[Act ^ Aa rn d] counts the number of Aa 
whose value is not smaller than Aa m d. As n — >• oo, it is clear that the estimated 
bootstrap p-value will tend to the ideal bootstrap p-value. Under the conventional 
significance level of 0.01, the multifractal phenomenon is statistically significant if 
and only if pi ^ 0.01. While pi > 0.01, the null hypothesis cannot be rejected. 

In a similar way, defining F = [f(a min )+f(a mBX )]/2, an analogous null hypothesis 
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can be described as follows: 
where the false probability is 



H*:F^F rndl (11) 



P2 ^ F > F ^ . (12) 
n 



Using the significance level of 0.01, the multifractal phenomenon is statistically 
significant if and only if j? ^ 0.01. 



3.3 Ensemble average 



In previous studies concerning multifractality in financial markets, the multifrac- 
tal analyses were carried out on individual stocks or index series. The index can 
be regarded as an ensemble average of the stock market in some sense. Here, 
we introduce a method to investigate the multifractality in an ensemble of many 
stocks, which is borrowed from the conce pt o f computing multifractal dimensions 
in Diffusion-Limited Aggregation 11461 I47U 48I. 14911 . Regarding the stock market as a 
stochastic process, a stock trading on the market can be considered as a realization 
of the stochastic process. We define quenched and annealed mass exponent TQ^(q) 
as follows 



m (XgO)) 



-r Q {q) Ins, 
-T A (q) Ins, 



(13) 
(14) 



where the angular brackets (•) is the ensemble average over all the chosen stocks. 
Intuitively, the annealed exponents are more sensitive to rare samples of the ensem- 
ble with unusual values of Xq( s )> while the quenched exponents are more charac- 



teristic of typical members of the ensemble [49J]. 



4 Results and discussions 



4.1 Multifractal analysis 



Two important stock indexes (SSEC and SZSC) and two stocks (Sinopec, 600028, 
in Shanghai Stock Exchange and China International Marine Containers, 000039, 
in Shenzhen Stock Exchange) are chosen as examples to show multifractality in 
single index or stock volatility series. Fig. [TJa) and (b) show the dependence of 
box size s for different values of q in log-log coordinates. Power 
laws with good quality are observed between Xq{ s ) l ^ q an d s. We also find that 
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the scaling range is wide and spans about three orders of magnitude. More interest- 



ingly, when s is small, there is a sudden jump on the Xq( s 



,1/(9-1) 



curve for q 



in Fig. QIb). This is not surprising for individual stocks since they may have time 
intervals within which the prices do not change. In this case, there is at least one 
box with very vanishing measure ji when the box size s is smaller than or equal 
to some critical value s c . When s > s c , the relative difference among measures is 
not large. Hence, for all negative q, we will observe such jumps. This explains the 
occurrence of a jump in individual stocks but not in indexes. Consequently, in the 
determination of the scaling ranges, we identify and exclude those jumps. 
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Fig. 1. (color online) Multifractal analysis, (a) Plots of Xq{ s ) 1 ^ q ~ l ^ > as a function of 
box size s for q = —3, 0, 3, and 5 in log-log coordinates for two indexes, (b) Plots of 
Xq(s) 1 ^ q ~~ 1 ^ as a function of box size s for q = —3,0,3, and 5 in log-log coordinates 
for two individual stocks, (c) Dependence of scaling exponents r(g) on the order q. (d) 
Multifractal spectra /(a) obtained by Legendre transform of r(g). 

There are at least three cases corresponding to "freezing" price in certain time pe- 
riod. The first case corresponds to the situation that the price has reached its daily 
price limit. According to the trading rules of SZSE and SHSE, there is a price fluc- 
tuation limit of ±10% for normal stocks and ±5% for specially treated (ST) stocks 
compared to the closing price of the last trading day. When the price reach its daily 
price limit with a huge number of shares waiting on the corresponding best bid or 
ask price, the price does not change for a long time, which might last till the closure 
of the market. The second case corresponds to those very liquid stocks for which 
both the buy and sell sides of the order book are very thick so that the price does 
not change frequently for most marketable orders. The third case corresponds to 
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those stocks with very low transaction activities. 

The scaling exponents r(q) are given by the slopes of the linear fits to In Xq( s ) with 
respect to In s for different values of q. We do this for all moments between -3 and 
5 with an increment of 0.2. Fig. [Qc) plots the dependence of the mass exponents 
r(g) on the moment order q. Fig. Old) presents the multifractal singularity spectra 
f(a) obtained through Legendre transformation of r(q) defined by Eq. ©. It is 
well-known that Aa = a mSLX — a min is an important parameter qualifying the width 
of the extracted multifractal spectrum. The larger is the Aa, the stronger is the 
multifractality. At the first glance, we find that Aa S SEC < Aa 6 ooo28 and Aa S zsc < 
Aaoooo39- This is in agreement with our common sense that the fluctuations of index 
are less volatile than that of individual stocks. 



To exhibit the analogue between financial market and fluid mechanics, the p-model, 
which is a simple cascade model of energy dissipation in fully developed turbu- 
lence 1150m . is applied to fit the mass exponent functions and it gives an excellent 
parametrization of the data. The theoretical mass function of the p- model can be 
expressed as follows 

\n[p q + (1 -p)i] 



r(q) 



In 2 



(15) 



The solid line shown in Fig. [TJc) is drawn according to the average of parame- 
ters obtained from fitting the four samples. We see that the agreement between the 
p-model and the stock data is remarkable for both positive and negati ve p arts, in- 
dicating the existence of information cascade process in stock market 115 ll. 15211 . We 
further extract the parameter p of the p-model for all the stocks under investigation. 
Fig.[2]illustrates that the empirical occurrence frequency g(p) as a function of p. We 
find that (p) = 0.40 ± 0.02. In contrast, fully developed turbulence gives p = 0.3 
& 
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Fig. 2. Empirical occurrence frequency g(p) of p. 
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4.2 Statistical tests 



We assess the statistical significance of the empirical multifractality in the sprit of 
bootstrapping tests. We reshuffle the time series and perform the same multifractal 
analysis. Fig. [3] compares the multifractal spectra of the raw time series and that 
of the 10 shuffled time series for the two chosen indexes and two securities. The 
thin lines are associated with the real data, while the thick lines are obtained from 
the shuffled data. An eye inspection already shows the deviation of the multifractal 
spectra of the real data f(a) and that of the shuffled data / m d(a m d)- We can infer 
that one of the most important causations for the existence of the multifractality is 
the long memory in the volatility series. 
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Fig. 3. Comparison between multifractal spectra extracted from real and shuffled data. The 
thin lines correspond to the real data, while the thick lines are for the shuffled data, (a) 
SSEC, (b) SZSC, (c) 000039, and (d) 600028. 

In the following, we shall give a more systemic statistical test on Chinese stock 
market. For all the chosen series, we shuffle the data 1000 times for each stock 
and reinvestigate the multifractality of the surrogate data. And the correspond- 
ing multifractal spectra are obtained. For each singularity spectrum, we calculate 
two characteristic quantities, Act and F = [f(a mill ) + /(a max )]/2. We find that 
Act > (Aa md ) and F < (F rnd ) for all the series, which implies that there are 
discrepancies between the multifractal spectra of the shuffled data and that of the 
real data. We find that p 1 = and p 2 = for the SSEC and pi = and p 2 = for 
the SZSC, which provides undoubtable evidence for the presence of multifractality 
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in the SSEC and SZSC data. Under the significance level of 1%, we find that the 
multifractal nature is significant for all stocks. 



4.3 Quenched and annealed average 



Fig. Ufa) illustrates the quenched and annealed exponents r 2 ^ against moment 
order q. For comparison, we also plot the mass exponents r(q) of SSEC and SZSC. 
The multifractal spectra obtained from Legendre transformation are presented in 
Fig. Kb). There are noticeable discrepancies between Aa^g and AassEC,szsc- 
Indeed, ensemble averaging allows us to capture the fluctuations among different 
realizations, which widens the singularity spectrum when compared with that from 
individual time series. This feature is well illustrated in Fig. 4(b). Moreover, one can 
see that Aa^ > Aag. The appearance of this inequality is linked to the sensitivity 
to rare values of the samples for annealed average. 
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Fig. 4. (color online) Quenched and annealed multifractal analysis, (a) Plots of mass expo- 
nents r(q) with respect to the moment order q. (b) Multifractal spectra. 



5 Concluding remarks 



We have performed detailed multifractal analysis on minutely volatilities of two in- 
dexes and 1 139 stocks in the Chinese stock markets based on the partition function 
approach. The minutely volatility is calculated as the sum of absolute returns in an 
interval of one minute with higher-frequency data. A measure fi is constructed as 
the normalized volatility. Hence, the measure /i is additive and conservational. This 
allows detailed multifractal analysis based on the partition function approach. We 
confirmed that \i is a multifractal measure. 

According to our analysis, the partition function Xq( s ) f° r ea ch security scales as 
a power law with regard to box size s for each order q. The function of scaling 
exponents r(q) for each security can be estimated. The nonlinearity of r(q) acts 
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as a hallmark of multifractality. The p-model in turbulence has been used to fit the 
r(q) function, resulting in a parameter p = 0.40 ± 0.02. Statistical tests based on 
the bootstrapping technique confirm the significance of the multifractal nature in 
the volatility time series of individual securities. 

An ensemble multifractal analysis was also carried out upon many stocks and the 
annealed and quenched mass exponent functions r(q) have been determined. We 
note that the ensemble averaging allows us to characterize the global multifractal 
properties of large ensembles of data (the market as a whole) in a compact form 
without dealing with details of individual realizations (individual stocks). In this 
sense, we can draw a conclusion that the Chinese stock market as a whole also 
exhibits multifractal behavior. 



As discussed in Section \T\ these correctly extracted multifractal characteristics 
might have potential usefulness in market prediction or risk management for in- 
dividual securities or the whole market. It is thus interesting to repeat the analysis 
in literature II32I. I341. 13811 based on our results. However, this is beyond the scope of 
the current work. 
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